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Abstract 



Using a recently developed mesoscopic theory of protein dielectrics, we have calculated the salt bridge 
energies, total residue electrostatic potential energies, and transfer energies into a low dielectric amyloid- 
like phase for 12 species and mutants of the prion protein. Salt bridges and self energies play key roles in 
stabilizing secondary and tertiary structural elements of the prion protein. The total electrostatic potential 
energy of each residue was found to be invariably stabilizing. Residues frequently found to be mutated in 
familial prion disease were among those with the largest electrostatic energies. The large barrier to charged 
group desolvation imposes regional constraints on involvement of the prion protein in an amyloid aggregate, 
resulting in an electrostatic amyloid recruitment profile that favours regions of sequence between alpha helix 
1 and beta strand 2, the middles of helices 2 and 3, and the region N-tcrminal to alpha helix 1. Wc found 
that the stabilization due to salt bridges is minimal among the proteins studied for disease-susceptible human 
mutants of prion protein. 
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Introduction 



Misfolded prion protein is the causative agent for a unique category of human and animal neurodegenerative 
diseases characterized by progressive dementia, ataxia, and death within months of onset (Prusiner 1998). 
These include Creutzfeldt- Jakob disease (CJD), fatal familial insomnia, and Gerstmann-Straussler-Scheinker 
syndrome in humans, bovine spongiform encephalopathy in cattle, scrapie in sheep, and chronic wasting 
disease in cervids. Unlike other infectious conditions that are transmitted by conventional microbes, the 
material responsible for propagation of prion diseases consists of an abnormally folded conformer of an 
endogenous protein, possibly in complex with host nucleic acids or sulfated glycans (Caughey et al. 2009). 
Soluble, natively- folded monomers of the prion protein (known as PrP^) may adopt an aggregated protease- 
resistant conformation known as PrP^'^that is capable of recruiting additional monomers of PrP^and inducing 
them to misfold in a process of template-directed conversion. This results in ordered multimers of prion 
protein that, when fractured, act as additional seeds to propagate the misfold through the reservoir of PrP*-^ 
present in brain. Although the conversion process may be initiated by an infectious inoculum of PrP^^, it 
may also arise spontaneously or due to mutations in the gene coding for PrP that predispose to misfolding. 

Structurally, PrP"-' is a glycophosphatidylinositol-anchored glycoprotein of 232 amino acids comprising 
an N-terminal unstructured domain and a C-terminal structured domain of 3 a-helices (hereafter referred 
to as al, a2, and a3 in order) and a short two-stranded antiparallel /3-sheet (made of strands (31 and (32), 
while PrP^'= has substantially enriched (3 content speculated to form a stacked ^-helix (Govaerts et al. 2004) 
or extended /3-sheet (Cobb et al. 2007) conformation in the amyloid fibril. 

At a molecular level PrP misfolding is a physico-chemical process, with the propensity to misfold deter- 
mined by the free energy difference between folded and misfolded states and the magnitude of the energy 
barrier separating them. As in any protein system, electrostatic effects make significant contributions to the 
energies of the various states and take two forms: salt bridge energy due to spatial proximity of charged 
groups within the native protein, and solvation/self energy due to field energy storage in the ambient protein 
and water dielectric media. A priori, it is expected that electrostatic effects generally favour the well-solvated 
monomeric PrP*-' over the more hydrophobic amyloid PrP^'^, since formation of PrP^*^ necessitates disruption 
of salt bridges in the native structure (although this may be compensated for by the formation of alternative 
salt bridges in PrP^'') and transfer of some charged groups into an environment of lower permittivity, both of 
which are energetically costly. However, these penalties on formation of PrP^'^ are counterbalanced by hydro- 
gen bonding, hydrophobic, and possibly entropic contributions that favour the amyloid form (Tsemekhman 
et al. 2007). Regional variation in the electrostatic transfer energy to water and amyloid may be useful in 
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predicting participation in the amyloid core of PrP"'^. Furthermore, several of the causative mutations for 
familial prion disease involve substitution of charged residues for uncharged residues (such as the D178N 
mutation responsible for fatal familial insomnia or familial CJD, depending on mutant allele polymorphism 
status at codon 129) or charge reversal of a residue (such as E200K, the most common mutation in classical 
familial CJD) (Kovacs et al. 2002), offering an indication of the importance of electrostatic effects in the 
misfolding process. More broadly, it has been found that changes in the charge state of a mutant protein 
compared to wild- type relate to its tendency to form aggregates (Chiti ct al. 2003), and the aggregation 
propensity of a polypeptide chain is inversely correlated with its net charge (Chiti et al. 2002); similarly, ag- 
gregation propensity is maximal at the protein iso-electric point where the net charge is zero (Schmittschmitt 
and Scholtz 2003). Intrinsically unstructured proteins tend to have a high net charge (Uversky, Gillespie 
and Fink 2000), which increases the electrostatic cost for the system to condense into the folded structure. 
Sequence correlations between charged groups may affect the kinetics of amyloid formation as well (Dima 
and Thirumalai 2004). 

The role of salt bridges in prion disease has been investigated previously by molecular dynamics simulation 
(MDS) and experimental studies of mutant protein. MDS of human PrP'-'has identified salt bridges that play 
a role in stabilization of the native structure (Zuegg and Grcady 1999). Other MDS studies of the R208H 
mutation, which disrupts a salt bridge with residues D144 and E146 of a-helix 1, have shown that it results 
in global changes to the backbone structure (Bamdad and Naderi-Manesh 2007). Experimentally, the E200K 
mutant of PrP'-'has been shown through calorimetry to be 4 kJ/mol less stable than wild type (Swietnicki et 
al. 1998). Mutation of two aspartates participating in al intra- helix salt bridges to neutral residues increases 
misfolding fourfold in cell-free conversion reactions under conditions favouring salt bridge formation (Speare 
et al. 2003). Interestingly, complete reversal of charges in al appears to inhibit conversion, possibly by 
preventing docking of PrP*-' and PrP^^^ (Speare et al. 2003). The pH dependence of charge interactions in 
PrP*"^ has also been investigated to identify those most sensitive to pH changes (Warwicker 1999); this is an 
important aspect of the problem because of the observed increased PrP*"^ misfolding rate at low pH. 

A unifying analysis of all PrP'-' salt bridges would be useful in understanding their role in structural 
stability. As well, to our knowledge solvation energy contributions to the misfolding process have not yet 
been investigated; they would offer a helpful perspective for probing the propensity of different regions of 
the prion protein to participate in the PrP^^ amyloid core. 

Direct extraction of salt bridge and solvation energies from molecular dynamics is complicated by the 
need to run long-length simulations that sample the equilibrium between states of interest, which can be 
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prohibitively slow for states that differ significantly in energy. An alternative approach is to use a continuum 
electrostatics description of the protein-water system, in which the response of surrounding material is 
modelled through solution of the Poisson-Boltzmann equation as a macroscopic dielectric that varies from a 
low value (usually 4) within the volume of the protein to 78 (the dielectric constant of bulk water) outside 
the protein. The downside of this method is that it ignores subtleties of the protein response to perturbing 
fields, such as cooperative internal reorganization. Using results from Kirkwood-Frolich theory (Ostcr and 
Kirkwood 1943; Frohlich 1949; Voges and Karshikoff 2000), we have recently developed a procedure to 
compute a spatially-varying dielectric function for a protein based on fluctuation statistics obtained from 
brief equilibrium MD simulations that capture much of the microscopic response of the protein at moderate 
computational cost (Guest et al. 2009). This provides a convenient tool to calculate solvation and salt 
bridge energies for all residues in a protein from a single simulation. In what follows we apply this method 
to determine the energies for all salt bridges in 12 molecular species of prion protein and the transfer energy 
for all residues in these proteins into a hypothetical protein amyloid core. 



Twelve structures of various species and mutants of PrP*-^were selected from the Protein Data Bank (PDB), 
including the species human IQLZ and IQLX (Zahn et al. 2000), cow IDXO (Lopez Garcia ct al. 2000), 
turtle 1U5L, frog IXUO, chicken 1U3M (Calzolai et al. 2005), mouse 1AG2 and IXYX (Rick ct al. 1996; 
Gossert et al. 2005), dog IXYK, pig IXYQ, cat IXYJ (Lysek et al. 2005), wallaby 2KFL (Christen et al. 
2009) and the human mutants D178N 2K1D (Mills et al. 2009) and E200K IFKC (Zhang et al. 2000). They 
were taken as starting points for 5ns all-atom molecular dynamics simulations using the CHARMM force field 
version C31B1 (Brooks et al. 1983) with explicit pure solvent water (no salt), periodic boundary conditions, 
particle mesh Ewald electrostatics, a timestep of 2fs, and a Lennard-Jones potential cutoff distance of 13. 5A. 
The basic residues (ARG and LYS) were protonated, while the acidic residues (HIS, ASP, and GLU) were 
deprotonated to reflect ionization conditions at pH 7. The system was first minimized for 200 time steps 
before starting the simulation. Snapshots of the simulations were taken every 2ps to build up an ensemble of 
equilibrium conformations for each protein. The dipole moments fj, of all residue side chains and backbones 
were calculated at each snapshot and used to obtain the correlation coefficients 
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for all pairs of Cartesian dipole components /Xj and /Xj, where the angle brackets denote an average over 
all snapshots (the thermal average). The matrix R of correlation coefficients was diagonalized to isolate 
the normal modes of dipole fluctuations, which describe the response of charged groups to perturbations 
around equilibrium. The R matrix for each protein was used to calculate the local dielectric map (Guest 
et al. 2009). See Figure 1 in the Supplementary Material for the dielectric map of human PrP. These 
dielectric maps were then taken as input for the Poisson-Boltzmann solver APBS (Baker ct al.2001) to solve 
the linearized Poisson-Boltzmann equation on an 97*^ mesh in 150 mM NaCl, again with periodic boundary 
conditions, to obtain the electrostatic energies required below. Atomic radii were assigned according to 
the CHARMM force field by the program PDB2PQR (Dofinsky et al. 2004). The often-used simplifying 
approximation of a constant internal protein dielectric constant of 4 and water dielectric constant of 78 was 
employed for comparison (Kumar and Nussinov 1999). 

Salt bridges in the set of proteins were identified by searching all pairs of charged atoms for those with 
charged groups within 12A of each other, whether the charges were alike or different. The energy of each 
salt bridge was determined by a mutation cycle designed to isolate the charge interaction energy from the 
energy in the surrounding dielectric milieu as shown in Figure 1. For charged groups A and B, their salt 
bridge energy Esb was taken to be a function of the energy of the protein system with both charges in place, 
Eab, with one or the other charge removed, Ea and Eb, and with both charges removed, Eq, as follows: 

Esb = Eq + Eab — {Ea + Eb) ■ (2) 

Here, Eq contains only the self energy of the part of the protein not including A and B (labelled P in Figure 
1), while Eab contains the self energies of A, B, and P as well as the pairwise interaction energies between 
A and B, A and P, and B and P. Ea contains the self energies of A and P and their interaction energies {Eb 
is analogous). Combining the terms as shown causes all the energies except the interaction energy of A and 
B to cancel. 

Another cycle, also shown in Figure 1, was used to determine the total contribution of each residue to 
the electrostatic energy of the protein, i?e;ec- For each side chain in the protein, the electrostatic energies of 

the side chain Egc and the protein lacking the side chain E^hoie-sc were calculated in isolation in the protein 
dielectric environment and subtracted from the electrostatic energy of the intact protein Eyj^oie- 

Eelec ~ Eyj^QiQ EgQ EyjfiQi^—gQ. (3) 
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The terms Eg^ and Eyj^oie-sc contain the self energies of the side chain and rest of the protein respectively, 
and Eyjhoie contains these self energies as well as the interaction energy of the side chain with the rest of 
the protein. Subtracting the terms as shown causes the self energies to cancel, leaving only the interaction 
energy between the side chain and the protein. This energy can be thought of as the electrostatic potential 
energy of a residue in the protein. 

To approximate the electrostatic energy of residue transfer into a hydrophobic, low dielectric environment 
like the core of a PrpSc amyloid, the energy E{eprpc) of a residue in the dielectric environment of PrP^ 
was compared to the energy Eiep^psc) of the residue in a homogeneous dielectric of e = 4, which describes 
the dielectric response in the interior of a bulk amyloid protein phase. Since the nature of monopole fields 
in the PrP^'= structure is unknown, interactions between charged residues are omitted from the calculation, 
so the transfer energy reflects only changes in the dielectric environment. For a given residue, the dielectric 
contribution to the transfer energy Etrans is: 

Etrans = E{epj.pSc) — E{eprpc). (4) 

Results 

Dynamics of Dipoles at Equilibrium 

The modes obtained by diagonalizing the correlation matrix R in Equation 1 generally involved several parts 
of the molecule; correlations were not limited to residues close in space or sequence. This is consistent 
with phonon transmission of perturbations at one site throughout the molecule by strong steric coupling 
effects through solid-like elastic moduli. The four largest-amplitude dipole modes for human PrP are shown 
in Figure 2. Dipole fluctuations were not qualitatively different between species, but different regions of 
the molecule exhibited characteristic motions. The two long alpha helices 2 and 3 exhibited primarily 
synchronous motion, with the helices rocking back and forth together as a unit. Nonetheless, some dipoles in 
the helices exhibited contrary motion. Alpha helix 1 did show some autonomy from the rest of the structure 
and tended to fluctuate as a group. 

Motion of the beta sheet is prominent in several of the modes. Two patterns stand out: a see-saw motion 
in which one strand tilts up as the other tilts down with both strands pivoting about the middle of the 
strands, and an in-out motion in which the outer beta strand (^1) and the N-terminal part of a2 move 
synchronously away from the inner beta strand (/32). The first motion is seen in modes 2 and 3 in Figure 
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2 above, while the second motion is seen in other lesser-ampUtude modes. This is compatible with NMR 
observations of the beta sheet, which show slow exchange between a range of conformations (Liu et al. 1999; 
Viles et al. 2001). In the NMR experiments, motion of the beta strands was observed on a time scale of 
microseconds, while these simulations only spanned nanoseconds, but both are indicative of some degree of 
conformational flexibility in the beta sheet. 

Salt bridge energies 

The PrP structures analysed contained a diverse set of salt bridges, ranging from moderately attractive to 
weakly repulsive. A complete list of salt bridges in all structures is presented in Table 1 of the Supplementary 
Material; salt bridges in the human structure are shown in Table 1 for both the single NMR structure IQLX 
and the ensembe of 20 NMR structures IQLZ. 

Structurally, the salt bridges can be divided into local and nonlocal by the proximity in sequence of the 
participating residues. Local salt bridges, like Aspl48 — Glul52 in al, Asp208 — Glu211 in a3, and Argl64 — 
Aspl67 between /32 and the following loop serve to stabilize secondary structural elements of the protein, 
while nonlocal salt bridges like Argl56 — Glul96, Argl64 — Aspl78, and Glul46 — Lys204 help to hold these 
elements together in the overall tertiary fold. Figure 3 shows the position of these nonlocal salt bridges in 
bovine PrP. 

Many of the salt bridges identified were near the protein surface, where the high degree of solvation 
attenuates their strength; the strongest salt bridges were those best sequestered from solvent, for this places 
them in a dielectric environment that increases electric field strength. The strongest salt bridge of all, 
between residues 206 and 210 of frog PrP, features a special "two-pronged" geometry that enables the 
amino group of Lys 210 to associate with both carboxyl oxygens on Asp 206. Interestingly, two strong but 
intermittant salt bridges arc present in human IQLZ between the C-tcrminal argininc and residue 167 in 
the /32 — a2 loop and residue 221 in aS. The substantial variation between members of the NMR ensemble 
at the C-terminus results in large motion of the arginine side chain, so that these salt bridges are only 
formed in a subset of conformers. Similarly, the ARG 164 - ASP 178 salt bridge that helps to anchor the 
beta sheet to a2 and a3 is not present in all members of the IQLZ NMR ensemble, although it is quite 
strong in the single IQLX structure. Although attractive salt bridges predominate, there were a number of 
repulsive salt bridges identified as seen on the left hand side of Figure 5A, especially in al and a3, which 
are crowded with several charged residues. As demonstrated in the following section, despite the presence of 
these destabilizing interactions, no residue experiences a net repulsive potential as these unfavourable salt 
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bridges are counterbalanced by the presence of other, stronger, favourable ones. The total energies due to 
all salt bridges in each molecular species studied are shown in Figure 5B. Of note is the much reduced total 
salt bridge energy in the two human mutants, E200K and D178N, compared to any other structure. The 
categorization of species as susceptible or resistant to prion disease is somewhat approximate, but comparison 
of total salt bridge energy and disease susceptibility by Kendall's tau gives a value of t = 0.45, implying 
that the order of species by salt bridge energy and disease susceptibility are significantly concordant (p = 
0.046). Overall, the effect of a heterogeneous dielectric was to moderate putatively strong salt bridges under 
the biphasic protein-water approximation for the dielectric function. 

The salt bridges listed in Table 1 of the Supplementary Material are those present at pH 7, but for human 
PrP an additional search was performed to identify salt bridges that would emerge at lower pH, since acidic 
conditions are known to drive PrP^*^ formation. Lower pH results in protonation of histidine residues to 
produce a positively charged species, which in human PrP enables the formation of three additional weakly 
attractive salt bridges (indicated by daggers in Table 1). While the dominant effect of lowering pH is to 
reprotonate acidic side chains, thus reducing electrostatic stability, this is partially compensated for by the 
formation of salt bridges involving histidine. 

Total residue electrostatic energies 

The salt bridge energies describe pairwise effects, but for mutational analysis it is more important to know 
the total contribution of each side chain to the stability of the protein. These energies approximate the 
electrostatic contribution to the energy change on mutation to a residue with a small nonpolar side chain 
like alanine. In practice, the side chain of each residue is removed from the protein. The total electrostatic 
energy of each residue in all prion proteins stiidied was less than or equal to 0, indicating a strong degree of 
evolutionary selection toward residues that benefit stability in the folded conformation. Through electrically 
neutral, or nearly so, proteins have their internal dipoles oriented so as to lower the potential energy of 
every residue. In human PrP, it is instructive to correlate the energies to known pathogenic mutations: 
the residue with the greatest overall stabilizing energy, Thrl83, is implicated in familial CJD by a T183A 
mutation (Kovacs et al. 2002); this mutation has also been shown to radically reduce measured stability 
by urea denaturation (Liemann and Glockshuber 1999). It is interesting to note that this residue, although 
not charged, is polar and more deeply buried in the hydrophobic core of the protein than any other charged 
residue, thereby enhancing the effect of dipolar attractions with its neighbors. Other residues that on 
mutation cause familial prion disease have especially high total electrostatic stabilizing energies, including 
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D178 and D202. Table 3 gives the 10 human side chains with the greatest total electrostatic energies. 
We might anticipate that mutation of other residues in Table 3 may enhance the probability of developing 
misfolding-related disease. 

Transfer to hydrophobic environment 

In forming the amyloid core of PrP^'^, some residues must undergo the migration to a region of low dielectric 
constant. For highly charged residues, this transfer energy is prohibitively high and may thereby exclude 
their participation in the amyloid core, while for nonpolar residues the small electrostatic transfer energy cost 
is overcome by favourable solvation entropy changes. By mapping the transfer energy of each residue into a 
region of low dielectric approximating PrP^'^ amyloid, it is possible to predict the likelihood of recruitment 
for various PrP regions into the amyloid core, without the aid of specific dipole-dipole correlations as might 
be present in the amyloid. Figure 6 shows the transfer energy from the PrP'-' dielectric to a homogeneous 
dielectric of 4 for various species of PrP^ . The transfer energy to an aqueous environment would show an 
inverse pattern. A 7 amino acid summing window is applied because sequence heterogeneity causes large 
variation between adjacent residues, and individual residues cannot enter the amyloid core without placing 
their neighbors in it as well. The transfer energies in Figure 6 are quite large, but including other terms 
in addition to the electrostatic energies considered here will reduce the magnitude of the total transfer free 
energy. 

There is considerable variation in the transfer energy along the sequence, with the lowest barrier to 
dielectric transfer for the region between al and 02, the middles of a2 and a3, and (31. Conversely, al, 
the loop between (32 and a2, and the loop between a2 and a3 show a formidable barrier to transfer. This 
overall pattern is well preserved among all PrP structures studied (see Supplementary Material Figure 
4). Immunological studies have defined (32 as a PrP^'^-spccific epitope (Paramithiotis et al. 2003), which 
presumably necessitates its surface exposure. In the human structure, (32 is located at the border between 
regions of low and high transfer energies, so it is possible that it is close in proximity to the amyloid core 
but protrudes sufficiently to be recognised by antibodies. 

The overall contour of the transfer energy functions is similar for all PrP structures studied, but there 
is some variation that correlates with known infectivity data. As seen in Figure 6, human and bovine share 
highly similar transfer energy profiles and are both susceptible to prion disease and interspecies transmission 
of disease, while non-mammalian turtle PrP that does not form PrP^'^ has a different profile, with a higher 
transfer energy barrier than cow or human over 4/5 of the sequence. PrP'-'from dog, a mammalian species 
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known to be resistant to prion infection (Polymenidou et al. 2008), is intermediate between the human 
and turtle profiles. The average transfer energies correlate with a species' resistance to disease (Figure 6 
legend). Some other species, including chicken, turtle, and wallaby, have a qualitatively different transfer 
energy function. 

Discussion 

Although electrostatic effects are only one contribution to the enigmatic PrP*" — > PrP^'^ conversion process, 
they offer clues to many of the central questions in prion biochemistry. The spatial variation of the dielectric 
is important, as neither charge separation distance, dielectric constant at the midpoint of the salt bridge, 
nor burial can predict salt bridge stability or total electrostatic energy (see Figures 2 and 3 of the Supple- 
mentary Material). As shown above, salt bridges play an important role in stabilizing both secondary and 
tertiary aspects of the PrP'-' structure. In fact, the total energy of all salt bridges in human PrP*-' is -60 
kJ/mol, almost twice the total stability of the protein as determined by calorimetry (Swietnicki et al. 1998; 
Liemann and Glockshuber 1999). Thus disruption of even a proportion of salt bridges in PrP'-' is sufficient 
to substantially destabilize the folded conformation, possibly accelerating or enabling the transition to PrP^'= 
in the right conditions. However, as has been observed elsewhere (Hendsch et al. 1996), the free energy 
change of salt bridge disruption may not be equal to the Coulombic energy of the salt bridge itself, due to 
the competing favourable reduction in desolvation energy. This will partially offset the change in stability 
from salt bridge disruption. Charge interactions may also participate in the poorly-understood association 
between the unstructured N-terminal domain and structured C-terminal domain of PrP^. As shown previ- 
ously, the C-terminal domain organization depends on the length of N-terminal tail present (Li et al. 2009), 
possibly through a collection of transient interactions below the detection threshold of NMR, resulting in an 
"avidity-enhanced" C-terminal structure. Charge complementarity between the N- and C-terminal domains 
provides one explanation for this phenomenon. For example, the very N-terminus of PrP*-' contains the highly 
positively charged region KKRPK from codons 25-29, while al contains the highly negatively charged region 
DYED from codons 144 147. If the N-terminal tail is free to explore a random walk around the C-terminal 
domain, electrostatic attractions are likely to bias this part of the tail toward residence near al, a region 
that is especially influenced by the length of tail present. The net attraction between the N-terminal tail 
and C-terminal structure may be insufficient to structure the tail but sufficient to collapse or condense the 
tail onto the surface of the structured domain, resulting in a kind of "molten shell." Further exploration of 
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this phenomenon, by molecular dynamics or other tools, may prove insightful. 

The importance of acidic conditions in the PrP conversion has been known for some time (Hornemann 
and Glockshuber 1998), and acidity exerts a large effect through modification of the protonation states 
of charged residues. At slightly acidic pH below the pKa of histidine (6.5), protonation of the histidine 
imidazole ring creates mildly stabilizing salt bridges with nearby residues. In the prion literature, histidines 
arc a subject of considerable attention for their ability to coordinate copper ions in the octapeptidc repeat 
region of the N-terminal domain (Aronoff-Spencer et al. 2000; Viles et al. 2008), but it seems that they also 
help to protect PrP^ from the stress of mildly acidic conditions. At much lower pH, however, protonation 
of glutamate and aspartate side chains ablates some of the stabilizing salt bridges shown in Table 1, which 
substantially reduces the energy barrier to rearrangement of PrP'-' components. For example, at pH 4.5, 
the pKa predictor program PROPKA (Bas et al. 2008) identifies glutamate residues 168, 200, 219, and 
221 as being significantly protonated, which will affect the stability of the protein and whose systematic 
investigation is a topic for future work. The influence of acidity on the monomeric PrP'-' structure has been 
extensively studied by molecular dynamics (Gu et al. 2003; Langella et al. 2004; DeMarco and Daggett 
2007), but perhaps the most noteworthy effect of acidity may not be in the resulting structural transition of 
the isolated PrP^ monomer but rather in lowering the barrier to induced reorganization in the presence of 
the templating species. 

Another natural question is the role that electrostatics play in the formation of the PrP^'^ amyloid. It 
may be argued that since the transfer energy profile in Figure 6 neglects the possibility of forming strong 
salt bridges in the low dielectric amyloid core of PrP^'^ it misrepresents the ability of these charges to stably 
occupy the amyloid. A counterexample may be constructed in a case of homogeneous dielectrics. The total 
energy change ^Etotal on bringing two opposite charges A and B both of charge q and radius Vion from a 
large distance apart in a medium with high dielectric esoiv hke water into close proximity tab in a region 
of low dielectric ep^ot to form a salt bridge is equal to the sum of the solvation energy changes Ai?^;^, and 
^■^foiv ^'^'-^ their pairwisc Coulomb energy /^Eab- Treating the charges as Born ions, in the limit where 
^prot/^soiv "C 1 (a valid assumption since generally egoiv = 78 and Cprot = 4, so Cprot/^soiv ~ 0.05), the total 
energy change to form the desolvated salt bridge is: 

AE,„,^, = AEi,^+AEg,^+AEAB = (J- - j^)--^ ^ 7^ i^- ' ^) ' 

ion / \^prot ^solv / ^prot' AB ^prot \'ion ' AB / 

This is always positive since tab is greater than rjon to satisfy the stereochemistry of the atoms. Thus 
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although salt bridges may partially mitigate the burial of charged residues, they cannot alter the fundamental 
unfavourability of the electrostatic component of this process. It is also possible that solvent-exposed salt 
bridges may form in the misfolded state outside or on the surface of the amyloid or oligomeric core, which 
could occur without the desolvation penalty described above. This would provide a mechanism to stabilize 
charged and polar parts of the protein in the misfolded form. Such salt bridges are likely to be relatively 
low in energy due to the high ambient dielectric environment, and it has been observed for amyloid-beta 
16-22 peptide that hydrophobic forces are more important than specific salt bridges in driving amyloid 
formation (Ma and Nussinov 2002). However as mentioned above, it has been shown that the net charge of 
a polypeptide chain incurs resistance to aggregation (Chiti et al. 2002; Schmittschmitt and Scholtz 2003), 
consistent with the notion of a higher overall energetic cost of transfer into a low dielectric medium for more 
highly charged polypeptides. In light of this, we believe the transfer energy profiles accurately convey this 
part of the obstacle to amyloid formation. 

Continuum electrostatics as a tool to examine protein behaviour has limitations, namely that it ignores 
the microscopic response of the system and thereby risks omitting subtle but important effects. However, by 
deriving the dielectric map from all-atom molecular dynamics simulations of the proteins of interest, we are 
able to substantially incorporate the microscopic response in our model and thereby improve the reliability 
of the energy estimates obtained. Previous theories could not reliably predict the effective dielectric constant 
inside a protein, so values typically between 4 and 10 have been used as initial guesses. Stronger salt bridges 
in the interior tend to be better predicted by an interior dielectric of 4, which would then overestimate the 
strength of the more abundant salt bridges on the protein surface. An interior dielectric of 10 best predicts the 
strength of the abundant surface salt bridges, but would then underestimate the strength of the buried interior 
salt bridges. The heterogeneous dielectric theory in (Guest et al. 2009) makes it unnecessary to guess at the 
value of the dielectric inside a protein and also indicates that no single value in the interior is satisfactory. 
Qiiantum effects due to electronic polarizability may be added to this approach as further refinement. The 
conformational variability in the ensemble of NMR structures for each PrP molecule also introduces an 
inherent uncertainty in the calculation of electrostatic energies, which we treated by averaging salt bridge 
energies over all NMR ensemble members. The molecular dynamics relaxation methods, often done in 
the absence of counterions, may introduce uncertainty as well. Electrostatic considerations are relevant to 
many aspects of the prion question, from PrP*-' dynamics and stability to PrP^'= amyloid organization and 
templating. We have presented an analysis of salt bridge, electrostatic, and hydrophobic transfer energies 
that provides a useful perspective for understanding the structural vulnerabilities of PrP'-'. 
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Esb ^AB Eq Eg 




Figure 1: Schematic of the approach for calculating salt bridge and total electrostatic energies. Circular 
arrows denote solvation/self energies; straight arrows denote interaction energies. Labels refer to quantities 
in Equations 2 and 3. 
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Figure 2: The 4 largest amplitude dipole correlation modes for human PrP. Regions of the same colour move 
in synchrony, while regions of different colours move in opposition. 




Figure 3: Stereo view of the three well-conserved nonlocal salt bridges as they are arranged in bovine PrP. 
The transparent surfaces show contours of equal dielectric (5 for blue; 70 for white) as determined from 
heterogeneous mesoscopic dielectric theory (Guest et al. 2009). The volume near the surface of the protein 
shows the greatest difference in dielectric on comparison of the homogeneous and heterogeneous dielectric 
fields. See also Figure 1 in the supplementary material for a surface plot of the dielectric permittivity as a 
function of position. 
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Frog 206 -210 (-29.7) Wallaby 1 56 - 202 (-1 3.3) Cat 156-202 (-1 1.8) Cat 147 - 1 51 (-1 1 .6) Chicken 159-215 (+5.2) 

Figure 4: The 4 most attractive salt bridges (in black lettering) and the 1 most repulsive (in red lettering) 
identified in the set of PrP structures. Numbers in parentheses are the salt bridge energies in kJ/mol for the 
structures indicated. Note that interactions with all surrounding residues were considered when assessing 
the total effect of each residue on overall stability, as given in Table 3 and Supplementary Table 2. 
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Table 1: Salt bridges in the human prion protein from IQLX (a single structure) and IQLZ (an ensemble of 
20 structures), f: Only present at low pH. The separation between charged groups is given by r. Esi,{l) is 
the salt bridge energy (an average for IQLZ) calculated with the heterogeneous dielectric theory, and £'sb(2) 
is the same energy calculated with a constant protein dielectric of 4. The standard deviation of salt bridge 
energy over the structures in the NMR ensemble containing the salt bridge is given as a fraction of the total 
salt bridge energy by The correlation coefficient between salt bridge energies from IQLX and IQLZ is 

0.82. 
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5 0-5 -10 -15 

Salt Bridge Energy (kJ/mol) 



^1 Resistant to prion disease 

Susceptible to infection but does not 
develop prion disease spontaneously 
I — I Susceptible to infection and develops 
— prion disease spontaneously 
^1 Familial prion disease mutations 




Figure 5: A. Histogram of average energies for all salt bridges identified. B. Total salt bridge energies in the 
molecular species studied. Error bars give the 95% confidence interval for the mean salt bridge energy from 
each ensemble of NMR structures. 
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Figure 6: Hydrophobic transfer energy for strands of 7 residues centred on a given residue index as calculated 
from Equation 4 for 4 species of PrP'-'. The numbers in the upper left hand corner give the average transfer 
energy over the whole protein. Shown below are the locations of secondary structural elements in the PrP^ 
sequence. 
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Table 2: The most attractive and repulsive salt bridges in the set of prion protein structures studied. The 
number of NMR conformcrs for each species in which the salt bridge is present is n. Esb{l) is the salt bridge 
energy (an average for IQLZ) calculated with the heterogeneous dielectric theory, and i?.,5(2) is the same 
energy calculated with a constant protein dielectric of 4. The standard deviation of salt bridge energy over 
the structures in the NMR ensemble containing the salt bridge is given as a fraction of the total salt bridge 
energy by 
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Table 3: Residues in human PrP (IQLX and IQLZ) with the greatest total eleetrostatic energy from Equation 
4. The last column gives the factor by which each residue's electrostatic energy exceeds the average for all 
residues in the protein (-33 kJ/mol). THR 183 has the lowest electrostatic energy in human PrP and is also 
a minimum on the hydrophobic transfer profile (see Figure 6), due presumably to its low dielectric local 
environment. VAL 210 appears in the list despite being a putatively nonpolar residue because it is located 
in a region of particularly low dielectric in the protein core, which increases the energy of its small side chain 
methyl group dipoles. /31-al: Between al and (31. al-(32: Between al and [32. Bold: Wild-type residues 
at the locations of known mutation sites in familial prion diseases. The correlation between transfer energies 
from IQLX and IQLZ is 0.87. With 90% confidence, the list is significantly enriched in pathologic mutations 
compared to random chance (p = 0.096). 



25 
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"Salt Bridges in the Prion Protein: 



Implications for Stability and Misfolding" 

o 




Figure 1: Local heterogeneous dielectric map of human PrP as calculated by mesoscopic dielectric theory (Guest et al., 
2009) . This is a slice taken through the geometric centre of the molecule in the X-Y plane. 



Table 1: Salt bridges in 13 structures of the prion protein with a 

repulsive or attractive energy greater than 3 kJ/mol. The number 
of NMR conformers in which the participating charged groups are 
within a distance of 12A is n. -Bsb(l) is the salt bridge energy cal- 
culated according to Equation 2 with the heterogeneous dielectric 
theory, and -Esb(2) is the same energy calculated with a constant 
protein dielectric of 4. The standard deviation of salt bridge en- 
ergy over the structures in the NMR ensemble containing the salt 
bridge is given by SEsb and as a fraction of the total salt bridge 
energy by 
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9.9 


0.66 


IQLZ 
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11 


-13.1 


0.9 


0.07 


-75.5 


5.7 


0.07 


IQLZ 
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12 


-13.7 


1.0 


0.07 


-76.6 


5.7 


0.07 
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Chicken 
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20 
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1.3 


0.25 
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7.3 


0.72 


1U3M 


Chicken 


155 159 


20 


-3.8 


0.7 


0.19 


-5.2 


2.3 


0.43 


1U3M 


Chicken 


171 185 


20 


-5.4 


1.6 


0.29 


-12.4 


6.7 


0.54 


1U3M 


Chicken 


159 215 


20 


5.2 


0.8 


0.15 


6.2 


2.7 


0.44 


1U3M 


Chicken 


163 219 


20 


-8.4 


4.3 


0.51 


-20.6 


15.9 


0.77 
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217 221 


20 


-5.7 


1.1 


0.20 


-15.3 


7.6 
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0.44 
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0.21 


-5.8 


2.3 
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1.42 
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0.55 
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-4.2 


2.8 


0.68 


-6.9 


7.3 
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0.43 
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0.46 
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-19.3 
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-5.3 
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11.9 


1.16 
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19 
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7.7 
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-21.4 


6.7 


0.32 
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15.9 


0.38 
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0.64 


-16.9 


13.6 
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224 


20 
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12.1 


1.05 
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-11.6 
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0.44 
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0.70 


IXYJ 


Cat 


146 


204 


20 
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0.32 
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0.52 
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0.51 


-9.3 


9.1 


0.97 


IXYK 
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147 


151 


20 


-8.4 


4.5 


0.54 


-22.5 


16.5 


0.73 


IXYK 
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148 


152 


20 


-3.1 


1.2 


0.39 


-5.3 


5.8 


1.09 
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0.60 
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1.25 
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-6.1 


0.8 
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0.38 


-5.6 
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1.11 


IXYK 
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20 
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0.7 


0.16 
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1.2 


0.26 


IXYK 
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20 


-5.3 


2.5 


0.48 


-11.5 


11.3 


0.98 


IXYQ 


Pig 


147 


151 


20 


-4.3 


1.3 


0.29 


-6.0 


4.1 


0.69 


IXYQ 


Pig 


148 


151 


20 


3.5 


0.3 


0.09 


4.1 


1.0 


0.24 


IXYQ 


Pig 


164 


178 


20 


-5.9 
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0.25 


-11.9 


7.0 


0.59 
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0.33 
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0.91 
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0.37 
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1.24 


IXYQ 


Pig 


220 


223 


19 


-3.8 


1.2 


0.33 
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0.67 
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0.25 
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0.38 


-9.2 
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0.44 
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1.18 
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0.28 
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0.61 
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3.9 
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0.46 
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0.17 
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1.2 


0.28 


2KFL 


Wallaby 


146 


208 


20 


-3.7 


0.3 


0.09 


-4.6 


0.9 


0.20 


2KFL 


Wallaby 


208 


211 


20 


-3.8 


0.9 


0.25 


-4.8 


5.0 


1.04 




' ' ■ • , I ' ' ' ' ' ' ' ' 
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Figure 2: A. Absolute salt bridge energies as a function of the reciprocal of the separation distance between the participating 
charged groups. B. Absolute salt bridge energies as a function of the reciprocal of the separation distance multiplied by 
the dielectric constant e at the midpoint of the line joining the charged groups. In a homogeneous dielectric environment, 
the energy of a salt bridge connecting charges qa and qb would be given by Est — j^fj^^TT^- Including 1/e in the fit in 
addition to 1/r improves the correlation from 0.63 to 0.76, but even taken together the quantities e and r do not reliably 
predict salt bridge energies. Thus a direct calculation of energies from the full heterogeneous dielectric map is necessary 
to obtain accurate results. 
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Table 2: Total electrostatic energies in kj/mol for every residue in a 
scl<;K;t.ion of 13 prion protein structures at pH 7. The energies are 
calculated according to Equation 3. 
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Figure 3: Total electrostatic energies of side chains in human PrP calculated from Equation 3 plotted as a function 
of solvent-accessible surface area (SASA), a statistic related to burial in the core of the protein. There is a weak but 
significant (p = 0.01) correlation between the energy and SASA, but SASA alone is a poor proxy for the electrostatic 
energy. 
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Figure 4: Hydrophobic transfer energies from Equation 4 for seven residue contiguous segments of 12 prion protein 
structures. The residue index gives the position of the middle residue in each segment. 
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